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ABSTRACT 

Radio-loud Active Galactic Nuclei at z ~ 2 — 4 are typically located in dense environments and their 
host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy 
evolution. Finding radio-loud quasars at the highest accessible redshifts (2 ~ 6) is important to study 
their properties and environments at even earlier cosmic time. They would also serve as background 
sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization 
in H I 21 cm absorption. Currently, only five radio-loud (R = /j,, 5 GHz//„ 44 ooA > 10) quasars are 
known at z ~ 6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan- 
STARRSl and radio FIRST surveys. The radio information allows identification of quasars missed 
by typical color-based selections. While we find no good 6.4 < 2 < 7.2 quasar candidates at the 
sensitivities of these surveys, we discover two new radio-loud quasars at z ~ 6. Furthermore, we 
identify two additional 2 ~ 6 radio-loud quasars which were not previously known to be radio-loud, 
nearly doubling the current z ~ 6 sample. We show the importance of having infrared photometry 
for 2 > 5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify 
the quasar J0203+0012 (z = 5.72), previously considered radio-loud, to be radio-quiet. Using the 
available data in the literature, we constrain the radio-loud fraction of quasars at z ~ 6, using the 
Kaplan-Meier estimator, to be 8.p3'2%- This result is consistent with there being no evolution of 
the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower 
redshifts. 

Subject headings: cosmology: observations - quasars: general 


1. INTRODUCTION 

The study of 2 ~ 6 quasars has shown the presence 
of almost complete Gunn Peterson absorption troughs 
in their spectra, corresponding to Lya absorption by the 
neutral Hydrogen in the intergalactic medium (IGM). 
This indicates a rapid increase in the neutral fraction of 
the IGM above z > 6, providing strong con straints on 
the e nd of the epoch of reionization (EoR) (|Fan et ah 
2006a |. The study of the IGM through Gunn-Peterson 
absorption has its limitations: quasar spectra suffer from 
saturated absorption at z > 6, and thus it becomes in- 
creas ingly difficult to study the IGM during the EoR 
(e.g., Furlanetto et al.|[2006 ). 

On the other hand, the 21 cm line (unlike the Lya 
line) does not saturate, allowing the study of the I GM 
even at large neutral fractions of Hydrogen (e.g., Carilli 
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et al.]|2004a). Therefore, the identification of radio-loud 


sources at the highest redshifts will be critical for current 
and future radio surveys. These objects will serve as 
background sources to study the intergalactic medium 
beyond the EoR through 21 cm absorption measurements 
(for example the Low-Frequency A rray and the Square 
Kilometer Array; e.g., see [Carilli et al. 2002). 

Typical high-redshift quasar searches are based on 
strict optical and near-infrared color criteria chosen to 
avoid the more numerous cool dwarfs which h ave si milar 
colors to high-redshift quasars (e.g., Fan et al(p 001). An 
alternative to find elusive quasars with optical colors in¬ 
distinguishable from stars is to require a radio detection. 
Most of the cool stars that could be confused with high- 
redshift quasars are n ot radio-bright at mJy sensitivities 
(Kimball et al. 2009); therefore, complementing an op¬ 
tical color-based selection with a bright radio detection 


reduces the contamination significantly (e.g., McGreer 
rtal.|009). - 


Currently, there are only three z > 5.5 quasars known 
with 1.4 GHz peak flux density > lmJy (J0836+0054, 
2 = 5.81, Fan et al. 2001; J1427+3312,' 2 = 6.12, Mc¬ 
Greer et al. 2006, Stern et al. 2007; and J1429+5447, 
2 = 6.18, Willott et al. 2010a). There are two other 
2 ~ 6 quasars in the literature classified as radio- 
loud, but with fainter radio emission (J2228+0110, 2 = 
5.95, Zeimann et al. 2011; and J0203+0012, z = 5.72, 
Wang et al. 2008). 

This small sample of currently identified radio-loud 
quasars at 2 ~ 6 has already provided important insights 
into Active Galactic Nuclei (AGN) and galaxy evolution, 
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emphasizing the importance of finding more of such ob¬ 
jects. For example, there is evidence that a radio-loud 
z ~ 6 quasar migh t be located in an overdensity of galax¬ 
ies (Zheng et al.|2006), similar to what has been found in 
radio -loud AGIN at lower redshifts (e.g.,|Wylezalek et al. 
20131. It has also been shown that another z ~ 6 radio- 
loud quasar resides in one of the most powerful known 
starbursts at z ~ 6 (O mont et al.||2013 ). 

At lower redshifts, it is well established that roughly 
10% - 20% of all quasars are radio-loud. It has been sug¬ 
gested that the radio-loud fraction (RLF) of quasars is 
a function of both optical luminosity and redshift (e.g ., 


Padovani 19931 La Franca et al. 1994||Hooper et al. 1995). 
If a differential evolution between radio-quiet and radio- 
loud quasars exists, it could indicate changes in prop¬ 
erties of black holes such as accretion modes, black hole 
masses, or spin (Rees et al.|1982 Wilson h Colbert|1995 


Laor 2000). This could provide insights into why some 
quasars have strong radio emission, while most have only 
weak radio emission. 

Some studies have found evidence of such evolution. In 


particular, Jiang et al. (2007) find that the RLF decreases 
strongly with increasing redshift at a given luminosity. 
For example, they find that the RLF at M 2 500 = — 26 
declines fro m 24% to 4% as r edshif t increases from 0.5 
to 3. Kratzer & Richards (2014) find a behavior in 
agreement with these hirelings m a similar redshift range 
(z ~ 0.5 — 2.5), but they also point out that the evolu¬ 
tion of the RLF closely tracks the apparent magnitudes, 
which suggests a possible bias in the results. However, 


little or no evidence of such evolution (e.g 

, | Goldschmidt 

et al.|1999t Stern et al. 2000[|Ivezic et al. 2 

002; Girasuolo 

eFaT 2003). 


In this paper we take advantage of the large area 


coverage and photometric information provided by the 
Faint I mages of the Radio Sky at Twenty cm survey 


(FIRST; Becker et al. 1995) and the Panoramic Survey 
Teles cope & Rapid Response S ystem 1 (Pan-STARRSl, 


PS1; Kaiser et al. 2002 2010) to search for radio-loud 


quasars at 5.5 > z A 7.2) We also revisit the issue of 
a possible evolution of the RLF of quasars with redshift 
by studying the RLF of quasars at the highest accessible 
redshifts, where an evolution (if existent) should be most 
evident. 

The paper is organized as follows. We briefly describe 
the catalogs used for this work in Section [5] The color se¬ 
lection procedures for 5.5 < z < 6.4 and z > 6.4 quasars 
with radio counterparts are presented in Section [3j In 
Section 4, we present our follow-up campaign and the dis¬ 
covery of two new radio-loud z > 5.5 quasars. The radio- 
loud definition used in this paper and details on how it is 
calculated are introduced in Section [5] In Section [ 6 j we 
investigate the radio-loud fraction of z >5.5 quasars by 
compiling radio information on all such quasars currently 
in the literature. This work identifies two additional 
high-redshift, radio-loud quasars that had not previously 
been noted to be radio-loud. 

We summarize our results in Section (TJ Magnitudes 
throughout the paper are given in the AB system. We 
employ a cosmology with Ho = 69 .3 km s -1 Mpc -1 , 
Hm = 0.29, and Ha = 0.71 (H inshaw et al.| 2013). 


2 .1. FIRST 

The FIRST survey was designed to observe the sky at 
20 cm (1.4 GHz) matching a region of the s ky mapped b y 


the Sloan Digital Sky Survey (SDSS; York et al. 2000) 
covering a total of about 10,600 square degrees. The 
survey contains more than 900, 000 unique sources, with 
positional accuracy to < 1". The catalog has a 5er de¬ 
tection threshold which typically corresponds to 1 mJy 
although there is a deeper equatorial region where the 
detection threshold is about 0.75 mJy. 

2 . 2 . Pan-STARRSl 

The PS1 3 - 7 T survey has mapped all the sky above decli¬ 
nation —30° over a period of ~ 3 years in five optical fil¬ 
ters c/pi, rp i, ipi, zpi, and ?/pi (Stubbs et al.|2010 Tonry 


et al. 2012). The PS1 catalog used in this work comes 
from the first internal release of the 37 r stacked catalog 
(PV1), whi ch is based o n the co-added PS1 exposures 


2. SURVEY DATA 


(see Metcalfe et al. 12013). This catalog includes data ob¬ 
tained primarily during the period 2010 May-2013 March 
and the stacked images consist on average of the co¬ 
addition of ~ 8 single images per filter. The 5cr me¬ 
dian limiting magnitudes of this catalog are gpi = 22.9, 
rpi = 22.8, ip\ = 22.6, zpi = 21.9, and yp\ = 20.9. PS1 
goes significantly deeper than SDSS in the i and z bands 
which together with the inclusion of a near-infrared y 
band allow it to identify new high-redshift quasars even 
in areas already covered by SDSS. 

3. CANDIDATE SELECTION 
3.1. The FIRST/Pan-STARRSl catalog 

In this paper we use the multiwavelength information 
and large area coverage of the FIRST and Pan-STARRSl 
surveys to find radio-loud high-redshift quasars. We 
cross match the FIRST catalog (13Jun05 version) and 
the PV1 PS1 stack catalog using a 2" matching ra¬ 
dius. This yields a catalog containing 334,290 objects. 
Given the similar astrometric accuracy of Pan-STARRSl 
and SDSS, we use the same matching radius u tilized b y 
the SDSS s pectroscopic target quasar selection (Richards 
et al.|2002). The peak of the SDSS-FIRST positional off¬ 
sets occurs at ~ 0 . 2 " a nd the fraction of false matc hes 
within 2" is about 0.1% (|Schneider et al.||2007H2010| see 
their Figure 6 ). Although this matching radius intro¬ 
duces a bias against quasars with extended radio mor¬ 
phologies (e.g., double-lobe quasars with out radio cores 
or lobe-donrinated quasars), Ivezic et al. (2002 1 show that 
less than 10% of SDSS-FIRST quasars have complex ra¬ 
dio morphologies. 

As redshift increases, the amount of neutral hydrogen 
in the universe also increases. At z > 6 the optically 
thick Lya forest absorbs most of the light coming from 
wavelengths A rest < 1216 A. This implies that objects at 
z ~ 7 (z ~ 6 ) are undetected or very faint in the z-band 
(i-band), showing a ‘drop’ in their spectra. They are 
thus called z-dropouts (i-dropouts). This ‘drop’ can be 
measured by their red z—y and i — z colors for z-dropouts 
and i-dropouts, respectively. 

Given that the radio detection requirement signifi¬ 
cantly decreases the amount of contaminants (especially 
cool dwarfs), we perform a much broader selection crite¬ 
ria in terms of colors and signal-to-noise ratio (S/N) in 
comparison with our z- and Adropout criteria presented 
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in Venemans et al. (2015) and Bahados et al. (2014), 
respectively, in particular for the z-dropouts in this pa- 
per, we allow for objec ts undetected in the zpi band to be 
0.1 mag bluer than in Venemans et al. (2015) and relax 
the S/N criteria in their q P1 and r P1 bands to be S/N < 5 
instead of S/N < 3 (see Section 3.2). For the i-dropouts 
in this p aper, we relax the selection limits in compari¬ 
son with Bahados et al. (2014). We allow the candidates 
to be 0.7 mag and 0.4 mag bluer in the ip± — zpi and 
rpi — zpi colors, respectively. For this selection we do 
not put any constrain in zpi — ypi color (see Section 
3.3). In this way we can detect 2 > 5.5 quasars which 
have similar colors to cool stars that are missed by typical 
color-based criteria. 


3.2. z-dropout catalog search (z> 6.4/ 

We require that more than 85% of the expected point- 
spread function (PSF)-weighted flux in the zpi and the 
ypi bands is located in valid pixels (i.e., that the PS1 
catalog entry has PSF_QF > 0.85). We require a S/N 
> 7 in the ypi band and exclude those measurements 
in the ypi band flagged as susp icious by the Image Pro¬ 
cessing Pipeline (IPP; [Magnier 2006] [2007) (see Table 6 
in Bahados et al. 2014). The catalog selection can be 
summarized as: 


S/N(2/ P1 )>7 (la) 

((S/N(z P1 ) > 3) AND (z P1 - ypi > 1.4)) OR 
((S/N(z P1 ) < 3) AND (z P1 , lim - yp\ > 1.3)) (lb) 
(S/N(i P1 ) < 5) OR (ipi - yp\ > 2.0) (lc) 
(S/N(r P1 ) < 5) (Id) 
(S/N(< 7 P1 ) < 5) (le) 

where zpiji m is the 3(7 zpi limiting magnitude. This 
selection yields 66 candidates. 

3.3. i-dropout catalog search (5.5 < z < 6.4/ 

Similar to the z-dropout catalog search, we require that 
more than 85% of the expected PSF-weighted flux in the 
ipi and the zpi bands is located in valid pixels. We 
require a S/N > 10 in the zpi band and exclude those 
measurements flagged as suspicious by the IPP in the 
zpi band. 

We do not put any constraint on the ypi band. This 
allows us to identify quasar candidates across a broad 
redshift range (5.5 < z < 6.4) and make better use of 
the Zpi band depth (which is deeper than the yp\ band). 
The yp\ information is used later on for the follow-up 
campaign. 

We can summarize the catalog selection criteria as fol¬ 
lows: 


S/N(zpi) > 10 (2a) 

((S/N(* P1 ) > 3) AND (i P1 - Zpi > 1.5)) OR 
((S/N(i P1 ) < 3) AND (t Plilim - z P1 > 1.0)) (2b) 
(S/N(rpi) < 3) OR (rpi - z P i > 1.8) (2c) 

(S/N( ffP1 ) < 3) OR (g P1 - z P1 > 1.8) (2d) 

where Jpyiim is the 3er ipi limiting magnitude. This 
query yields 71 candidates. 

3.4. Visual Inspection 


T he n umber of candidates obtained from Sections 13.21 
and 3.3 are small enough to visually inspect all of them. 
We use the latest PS1 images available for the visual in¬ 
spection, which are usually deeper than the images used 
to generate the PS1 PV1 catalog. We also perform forced 
photometry on them to corroborate the catalog colors (as 
described in Banados et al. 2014), especially when PV1 
only reports limiting magnitudes. Thus, we visually in¬ 
spect all the PS1 stacked, FIRST, and zpi and yp\ PS1 
single epoch images (where the S/N is expected to be the 
highest) for every candidate. The most common cases 
eliminated by visually inspection are PS1 artifacts due 
to some bad single epoch images, objects that lacked in¬ 
formation in the PV1 catalog, and objects with evident 
extended morphology in the optical images. Based on 
the visual inspection, we assign priorities for the follow¬ 
up. Low-priority candidates are the ones whose PS1 de¬ 
tections look questionable, in the limit of our S/N cut, 
and/or objects with extended radio morphology which 
produces slightly positional offsets (> 1") between the 
optical and radio sources. 


3.4.1. z-dropouts 

Almost all the candidates can be ruled out by their PS1 
stack images and/or single epoch images. There are only 
two objects we cannot completely rule out although there 
are some lines of evidence pointing us to believe they are 
unlikely to be quasars. For both PSO J141.7159+59.5142 
and PSO J172.3556+18.7734, the ypi detection looks 
questionable, in the limit of our S/N cut: 7.5 and 7.0, 
respectively. Their ypi catalog aperture and PSF mag¬ 
nitudes differ by 0.3 and 0.28 mag which could indicate 
they are extended sources but it is hard to tell at this low 
S/N. The optical and radio positional offsets are some¬ 
what larger than for most of the candidates: 0.6" and 
1.5". All of this combined makes them low priority can¬ 
didates. The PS1 and FIRST information for these 
sources is listed in Table [U 

3.4.2. i-dropouts 


One of the candidates we selected is the known radio- 
loud quasar ,10836+0054 at z = 5.82 (Fan et al.||2001|). 
Its images look good and we would have followed it up. 
There are two low-redshift quasars that could have been 
selected for follow-up, J0927+020 3 a quasar with a bright 


Schneider et al. 


201 0 ) 


Ha line at ~ 9200A (z = 0.39; 
and J0943+5417, an Fell l ow-ionization broad absorp- 
tion line quasar (z = 2.22; Urrutia et al. 2009). After 
the visual inspection and literature search, there are 10 
remaining candidates, out of which 9 are high priority 
candidates. Their PS1 and FIRST photometry are listed 
in Table H] 


4. FOLLOW-UP 
4.1. Imaging 


We use a variety of telescopes and instruments to con¬ 
firm the optical colors and to obtain near-infrared pho¬ 
tometry of our candidates, thereby allowing efficient re¬ 
moval of interl opers. __ 

GROND (Greiner et al. 2008) at the 2.2m telescope 
in La Silla was used to take simultaneous images in the 
filters grizJHK during 2014 January 24 - February 5. 
Typical on-source exposure times were 1440 s in the near- 
infrared and 1380 s in the optical. The ESO Faint Object 
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Spect rograph and Camera 2 (EFOSC2; Buzzoni et al. 
19841 and the infrare d spectrograph and ima ging camera 
Son of ISAAC (SofI; |Moorwood et al.||1998|) at the ESO 
New Technology Telescope (JNTT) were used to perform 
imaging in the Intt (*#705), ^ntt (z# 623), and Jntt 
bands during 2014 March 2-6 with on-source exposure 
times of 600 s in the Intt and Zntt bands and 300 s 
in the Jntt band. The data reduction consisted of bias 
subtraction, flat fielding, sky subtraction, image align¬ 
ment, and stack ing. The photometri c zeropoints were 


determined as in Banados et al. (2014 9 and their errors 


are included in the magnitudes reported in this work. 
All of our high-priority candidates were photometrically 
followed up except for one which we directly observed 
spectroscopically (see next Section). Two low-priority z- 
dropouts and one low-priority i-dropout are still awaiting 
follow-up. Table[2]shows the follow-up photometry of our 
candidates. 

4.2. Spectroscopy 

We have taken spectra of four high-priority objects 
that were not rejected by the follow-up photometry. We 
processed the data using standard techniques, including 
bias subtraction, flat fielding, sky subtraction, combi¬ 
nation of individual frames, wavelength calibration, and 
spectrum extraction. T he spectra were flux calibrated 
using standar d stars from Massey & Gronwall (1990) and 


Harnuy et al. 
The candle 


(1992 1994). 


ate PSO JT14.6345+25.6724 was observed 


with the FOcal Re d ucer /l ow dispers ion Spectrograph 2 
(FORS2; lAppenzeller & Rupprecht |l992) at the Very 
Large Telescope (VLT) on 2014 April 26 with an expo¬ 
sure time of 1497 s. The spectrum shows no clear break 
in the continuum and we classify this object as a radio 
galaxy at z = 1.17 by the identification of the narrow 
[O n] A3728 emission line. 

We obtained an optical spectrum of the candidate 
PSO J354.6110+04.9453 using the Double Spectrograph 
on the 5 m Hale telescope at Palomar Observatory 
(DBSP) on 2014 July 21 for a total integration time of 
1800 s. The object has a red spectrum lacking the clear 
(Lya) break which is a typical signature of high-redshift 
quasars. The spectrum does not clearly identify lines to 
determine a redshift. We believe this object is most likely 
a radio galaxy, but it is definitely not a z > 5.5 quasar. 
The other two spectroscopically followed up objects - 
PSO J055.4244-00.8035 and PSO J135.3860+16.2518 - 
were confirmed to be high-redshift quasars. These two 
newly discovered quasars would not have been selected as 
candidates by the o ptical selection criteria presented in 
Bafiados et al.|([20l4|) (although PSO J055.4244-00.8035 
has only a lower limit of Jpi — Zpi >1.3 and it might 
have been selected if deeper ipi data was available). The 
observations for these quasars are outlined in more detail 
below. 


time of 2400 s. These observations were carried out in 
~ 1" seeing using the 175 wide longslit. This spectrum 
shows a sharp Lya break indicating that the object is 
unambiguously a quasar at z > 5.5 but the S/N does not 
allow us to determine an accurate redshift. We took a 
second spectrum with FORS2 at the VLT on 2014 Au¬ 
gust 4; the seeing was l'.'l and it was observed for 1467 s. 
This spectrum is shown in Figure [l] and there are no ob¬ 
vious lines to fit and use to derive a redshift. There 
is, however, a tentative Si IV+ O iv] line which falls in 
a region with considerable sky emission and telluric ab¬ 
sorption making it not reliable for redshift estimation. 
We estimate the redshift instead by comparing the ob¬ 
served quasar spectrum with the com posit e SDSS z ~ 6 
quasar spectrum from Fan et al. (2006b). We assume 
the redshift that minimizes the x i between the observed 
spectrum and the template (the wavelength range where 
the minimization is performed is A res t = 1240 - 1450 A). 
The estimated redshift is 2 = 5.68. However, because of 
the lack of strong features in the spectrum, the y 2 distri¬ 
bution is relatively flat around the minim um and thus a 
range of redshifts is acceptable. We follow Banados et al. 
(20141 and assume a redshift uncertainty of 0.05. 


4.2.2. PSO J135.3860+16.2518 (z = 5.63 ± 0.05) 

The discovery spectrum was taken with EFOSC2 at the 
NTT on 2014 March 3. The observations were carried out 
with the Gr#16 grism, 1"5 slit width, 173 seeing, and a 
total exposure time of 3600 s. The spectrum is very noisy 
but it resembles the shape of a high-redshift quasar with 
a tentative Lya line at ~ 8100 A. In order to increase 
the S/N and confirm the quasar redshift, we took two 
additional spectra and combined them. One spectrum 
was taken on 2014 April 5 wit h the Multi-O bject Double 
Spectrograph (MODS; |Pogge et al.||2010|) at the Large 
Binocular Telescope (LBT). The observations were car¬ 
ried out under suboptimal weather conditions with the 
172 wide longslit for a total exposure time of 2400 s. The 
second spectrum was taken with FORS2 at the VLT on 
2014 April 26 with a total exposure time of 1467 s. The 
observing conditions were excellent with 0755 seeing and 
we used the 173 width longslit. The combined MODS- 
FORS2 spectrum is shown in Figure [T| The estimated 
redshift by mat ching to the z ~ 6 quasar composite spec¬ 
trum from Fan et al. (2006b) is z = 5.63 ± 0.05. The 
redshift estimate for PSO J 135.3860+16.2518 is quite un¬ 
certain as represented by its error bar. There might 
be a tentative Siiv+Oiv] line which would place this 
quasar at a slightly higher redshift. However, this line 
falls in the same region as the tentative line in the quasar 
PSO J055.4244-00.8035, which is not a reliable region for 
redshift determination. A higher S/N optical spectrum 
and/or a near-infrared spectrum would be beneficial to 
obtain a more accurate redshift. 


4.2.1. PSO J055.42U~00.8035 (z = 5.68 ± 0.05/ 

The discovery spectrum was taken on 2014 February 
22 using the DBSP spectrograph with a total exposure 

9 Color conversions missing in |Banados et al.| |2014|| : 
ffGROND = Spi + 0.332 X (g P1 - V P1 ) + O.Ubb; 
rGROND — r Pi + 0.044 x (r P i — ipi) — 0.001; and Jntt and 
JCqrOND are calibrated against 2MASS. 


5. RADIO-LOUDNESS 


A clear consensus on a boundary between radio-loud 
and radio-quiet quasars has been difficult to achieve and 
there are several radio-loudness criteria in the litera¬ 
ture (for a comparison of the different criteria see Hao| 


et al. 2014|. We adopt the most widespread definition 
in the literature. This is the radio/optical flux density 
ratio, R = //,5GHz// p 440 oA (iKellermann et al.||1989h, 
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Figure 1. Spectra of the two radio-loud quasars discovered in this 
paper. The gray solid lines show the lcr error in the spectra. The 
blue dashed line is the SDSS z ~ 6 composite quasar spectrum 
from |Fan et al.| (|2006b|) at the redshift of the quasars for compari¬ 
son. Vertical dotted lines indicate the observed wavelengths of key 
spectral lines, as given in the top panel. 


where /„, 5 GHz is the 5 GHz radio rest-frame flux density, 
fv 4400 A i s 4400 A optical rest-frame flux density, and 
a quasar is considered radio loud if R > 10. 


5.1. The Radio Emission 

The rest-franre 5 GHz radio flux density is obtained 
from the observed peak flux density at 1.4 GHz. We use 
the peak flux density since most of the z ~ 6 quasars 
appear to be unr esolved on t he radio maps ( e.g., Wang 
et al. 12007 2011[ but see also Cao et al.|2014| where they 
claim that there may be extended structures around the 
radio-loud quasar J2228+0110 on arcsecond scales). 

We assume a power-law (/„ ~ v a ) radio spectral en¬ 
ergy distribution, adopting a typical radio spectral index 
a« = —0.75 as used in other high-reds hift studies (e.g., 
Wang et al. 2007 Momjian et al.||2014|. This ind ex ap¬ 


pears to be appropriate as|hrey et al.|(|2005[ |2008[ 
find that at least three z ~ 6 radio-loud quasars s 
steep radio spectrum. 


2011 ) 


how a 


5.2. The Optical Emission 

The optical spectral indices of quasars span a fairly 
large range (at least — 1 < a v < 0). When a direct mea¬ 
surement of the optical rest-frame flux density at 4400 A 
is not possible, this is typically extrapolated from the AB 
magnitude at rest frame 1450 A (TO 1450 ) assuming 
erage optical spectral index of a v = —0.5 (e.g., 


an av- 
Wang 


et al.|2007). These fairly large extrapolations could lead 


to dramatic errors if the studied quasars are not average 
quasars. At z ~ 6 — 7 we can take advantage of infrared 
space missions such a s Spitzer and the Wide-Field In¬ 
frared Survey (WISE; Wright et al. 2010) to reduce the 
extrapolation error by using ~ 3 /rrn (A res t = 4286A at 
z = 6 ) photometry. 

Following previous works, we also assume an optical 
spectral index of a v = —0.5 but we estimate f v 4 400 a 
from IRAC ( Fazio et al.||2004 ) 3.6/im ( 63.6 pm) observa¬ 
tions when available (A res t,eff = 5071A at z = 6 ). Other¬ 
wise, f v 4400 a is calculated from the WISE W\ magnitude 


(A r est, eff = 4811A at z = 6 ) reported in the ALLWISE 
Source catalog or reject table (with S/N > 2.5). Table [3] 
lists S 3 . 6 Wi, and 7714450 for all 2 > 5.5 quasars witn 
published measurements at 1.4 GHz in the literature. 
There are five quasars without IRAC or WISE mea¬ 
surements. For these objects, f v 440 oA i s estimated from 
7774450 - To estimate the error we determine /„ 440 oA from 


7774450 f° r the quasars having IRAC or WISE (~ 3/im) 
measurements. Then, for each object we compute the 

ratio (/j /j4 4ooa( 3M m ) - A,44ooA( m i45o))/A,44ooA( 3 M™). 

This results in a symmetric distribution, centered on 
zero, and with a standard deviation of 0.4. Finally, we 
take the absolute value of the previous distribution and 
assume its median as a representative error for /„ 44 0 oA 
derived from 7774450 . This corresponds to a relative error 
of 0.30 (i.e., analogous to assuming a measurement of 
fv 4400 A w hh S/N = 3.3). These errors must be taken 
with caution since they are just representative uncertain¬ 
ties and there could be objects with considerably larger 
errors, as exemplified in Section | 6 T] 


6. RESULTS 

We determine the radio-loudness parameter R and the 
optical and radio luminosities L 4400 A an<3 Tsghz (L = 
vL v ) for all the z > 5.5 quasars having 1.4 GHz data pub¬ 
lished in the literature. These parameters are included 
in Table [3] Figure [2] shows the rest-frame 5 GHz radio 
luminosity versus rest-frame 4400 A optical luminosity. 
Eight quasars are classified as radio-loud (R > 10) , in¬ 
cluding the two quasars discovered in Section |4.2| and 
two additional possib le ra dio-loud quasars which will be 
introduced in Section [672] There are 33 objects robustly 
classified as radio-quiet quasars (R < 10). Two quasars 
need deeper radio data to classify them unambiguously 
(with radio loudness upper limit > 10). In this paper we 
do not find any radio-loud quasar at z ~ 7 in an area 
of about 10, 600 square degrees of sky to the sensitivities 
of the FIRST and Pan-STARRSl surveys (~ lmJy and 
7cr-limiting magnitude 7/pi ~ 20.5, respectively). Con¬ 
clusions on the RLF at z > 6.5 are not possible at this 


at z > 6.5 (Mortlock et al. 

2011 

Venemans et al. 

2013 

2U15). From these seven quasars, only one has dedicated 


is in th e FI RST footpr int but It i s~a non-detection (see 
Section |6.3| |Venemans et al.|2015|) . The RLF at 2 6 is 

i • - f 1-1 1-7-i r-s - I a a I- - 1 - 


discussed in Section 16.31 


6.1. J0203+0012: A Radio-Loud Quasar? 
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The quasar J0203+0012 at z = 5.72 jwas classified 


as radio-loud quasar by Wang et al. (2008). They de¬ 
rived the rest-frame 4400 A flux density from the 1450 A 
magnitude. As discussed previously, these large extrap¬ 
olations can carry large uncertainties and can be criti¬ 
cal for the classification and derived parameters for spe¬ 
cific objects. This is the case for J0203+0012, a broad- 


absorption line quasar (Mortlock et al. 2009) whose 
mi 45 o could have been underestimated, resulting in a low 
optical luminosity (see also i ts spectr al energy distribu¬ 


tion in Figure 14 of Leipski et al. 2014). This quasar has a 
radio-loudness parameter K = 4.3±0.5 or R = 12.1 ±3.2 
depending if the rest-frame 4400 A flux density is extrap¬ 
olated from the IRAC 3.6 urn photometry or from 7774450 , 
respectively (see Figure [2]). It is clear that the differ¬ 


ence is dramatic and thatby using the 7711450 proxy this 
quasar would be classified (barely) as radio-loud. We ar¬ 
gue that the value of L 4400 A obtained from the observed 
3.6/rm photometry is more reliable for z > 5.5 quasars 
since it relies less on extrapolation (for this particular 
case, the extrapolation is less by approximately a factor 



6.2. Pushing the FIRST Detection Threshold 

The FIRST survey has a typical source detection 
threshold of 1 mJybeam -1 which assures that the cat¬ 
alog has reliable entries with typical S/N greater than 
five. There are 30 known z > 5.5 quasars that are not 
in the FIRS T catalog o_r_in the Stripe 82 VLA Survey 


catalog (Hodge et al. 2011) but that are located in the 


are 

Mahabal et al. (2005): 

Cool et al. 

2006): Jiang et al. 

(20 

19 

i; 

); 

Wiilott et al. (2009 

201031 

): McGreer et al. 

(2013 

Banados et al. 

(2014 

); Venemans et al. (2015|; 


and Warren et al. (in prep.). 

We checked for a radio detection beyond the FIRST 
catalog threshold as follows: we obtained the 1.4 GHz 
FIRST images for all 30 quasars and checked for radio 
emission within 3" of the optical quasar position with 
a S/N> 3. We find that the quasars J1609±3041 at 
z = 6.14 (Warren et al. Jin prep) and J2053±0047 at 
z = 5.92 (Jiang et al. 2009) have tentative 1.4GHz de¬ 


tections at S/JNI of 3.5 and 3, respectively. Their radio 
postage stamps are shown in Figure [3] The two quasars 
have optical-to-radio positional differences less than 1"8 
(1 pixel). In order to quantify the probability of finding a 
spurious association with a 3<r fluctuation given our sam¬ 
ple of 30 quasars we performed the following steps. We 
placed 100 random positions in each 1 arcmin 2 FIRST 
image centered on a quasar and measured the maximum 
peak flux within T/8. We removed points falling within 
1'.'8 from an optical source in the PS1 catalog. There are 
no radio sources in the FIRST catalog for any of these 
1 arcmin 2 fields. We computed the fraction of measure¬ 
ments with S/N> 3. We repeated this procedure 100 
times and the fraction of measurements with S/N > 3 
was always < 1%. The full distribution is centered on 


0.5% with a standard deviation of 0.1%. Therefore, in 
our sample of 30 quasars, the expected number of spuri¬ 
ous ~ 3<r associations within 1"8 is 0.15 and being con¬ 
servative less than 0.3. This analysis suggests that these 
identifications are unlikely to be spurious. If these de¬ 
tections are real, these quasars are classified as radio-loud 
with R = 28.3±8.6 and R = 44.1 ±18.7 (see Table [3]). In 
Figure[2j these new tentative radio detections are marked 
as red downward triangles. 

We make a mean stack of the 28 remaining quasars in 
the FIRST footprint that have S/N < 3 in their individ¬ 
ual images. We find no detection in the stacked image 
with an upper limit of /„ = 84 /iJy (see Figure [ 2 ]). 


6.3. Constraining the Radio-Loud Fraction of Quasars 

at z ~ 6 

Considering all the quasars in Table [3j there are eight 
known radio-loud quasars at z ~ 6, 32 radio-quiet (ex¬ 
cluding J1120±0641 at z = 7.08), and two ambiguous. 
There is one additional quasar that is robustly classified 
by a non-dete ction in FIRST as radio-quiet: J0148±060 0 


at z = 5.98 (Banados et al. 
This radio-quiet quasar has 


2014 Becker et al. |2015| ). 
■°g A 4400 A (A©) = I3.04± 
0.2, log± 5 Ghz A©) < 8.7, and R < 5.6 (see Figure [2]). 
We can provide a rough estimation of the radio-loud frac¬ 
tion of quasars at z ~ 6 of RTF = 8/(8±(34±l)) ~ 19%. 
In this statistics we considered the two ambiguous cases 
as radio-quiet. This is a relatively large fraction; how¬ 
ever, these quasars were selected by several methods 
which could potentially bias the results. As we have in¬ 
cluded radio-loud quasars that could not have been dis¬ 
covered based on their optical/near-infrared properties 
alone, the actual fraction of radio-loud quasars is over¬ 
estimated. Therefore, this value has to be taken only as 
an upper limit. 

In order to set a lower limit in the RLF at z ~ 6 , we 
consider quasars that were selected based on their optical 
properties (i.e., we exclude the two quasars discovered in 
this paper, J2228±0110 which was discovered by its ra¬ 
dio emission by Zeimann et al. 2011, and J1427±3312 
which was discovered by its radio emission by McGreer 
et al. 2001]^]). We also exclude quasars at z > 6.5, i.e., 
J1120±0641 at z = 7.08 and PSO J036.5078±03.0498 


at z = 6.527. The latter quasar was discovered in Ven- 
emans et al.| (|2015|). It is not detected in FIRST and 
has log ± 44oo a( l 0 ) = 12.87±0.03, logL 5 G hz(A©) < 8 . 8 , 
and R < 10.7 (see Figure[2]). Considering all the FIRST 
non-detections from the previous section as radio-quiet, 
we find a lower limit of RLF = 4/(4 ± (34 ± 27)) ~ 6 %. 
This is a lower limit because there is still the possibility 
that a fraction of the FIRST non-detections are radio- 
loud (see Figure [2j. Therefore, in this case we are poten¬ 
tially underestimating the number of radio-loud quasars. 

Additionally, in order to fully use the information pro¬ 
vided by both the radio detections and upper limits, we 
estimate t he radio-loud fraction u sing the Kaplan-Meier 
estimator (Kaplan & Meier 1958). The RLF estimated 
with this method, after excluding quasars at z > 6.5 and 
quasars selected by their radio emission, is 8 .ljjg. 2 %. 


10 We n ote however , that J1427+3312 was independently dis¬ 
covered by|Stern et al. (2007J without using the radio information, 
but using a mid-infrared selection instead. 
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6.4. What Changes with an Alternative Radio-Loudness 
Definition? 

Another common definition for radio-loudness in the 
literature (besides our adopted criteria in Section [5]), is a 
simple cut on rest-frame radio luminosity. This criteria 


radio luminosities are not correlated ( 

Peacock et al. 

1986 

Miller et al. 

1990 

Ivezic et al. 

2002 

Appendix G 

). We 


fixed radio luminosity as a boundary between radio-quiet 
and radio-loud objects . We use the alternative criteria 
adopted by | Jiang et al.| (|2007|), where a radio-loud quasar 
is defined with a luminosity density at rest-frame 5 GHz, 
Ai/, 5 Ghz > 10 32 ' 5 ergss -1 Hz -1 . This is equivalent to re- 
quiringlogL 5 Qhz(^©) > 8.61 (see the horizonal line in 
Figure[2]). 

One caveat with this definition is that J2228+0110 
is just below the radio-loud cut although is still con¬ 
sistent with being radio-loud within the uncertainties: 
logL 5 Ghz(-b©) = 8.59 ± 0.08. Note that the quasar 
J0203+0012, discussed in Section |6.1[ is also classified 
as radio-quiet by this definition. 

We estimate the radio-loud fraction using the Kaplan- 
Meier estimator, following the approach of the previ¬ 
ous section. The estimated RLF with this definition is 
6.6j^K %. T his result agrees with the one obtained in 
Section |6.3[ and they are both consistent with no strong 
evolution of the radio-loud fraction of quasars with red- 
shift. 


7. SUMMARY 

We perform a search for high-redshift, radio-loud 
quasars (i- and 2 -dropouts) by combining radio and opti¬ 
cal observations from the FIRST and Pan-STARRSl sur¬ 
veys. The multiwavelength information of these surveys 
allows the identification of quasars with optical colors 
similar to the more numerous cool dwarfs and therefore 


missed by typical color selection used by high-redshift 

quasar surveys (e.g., Fan et al. 2006c 

Banados et al. 

2014). We do not find good quasar candi 
(^-dropouts). We discover two of th 
quasars at 2 > 5.6: PSO J055.4244-00. 
with a radio-loudness parameter R = 1 
PSO J135.3860+16.2518 (2 = 5.63) with 
These two quasars are at the low-redsh 
dropout selection technique (5.5 < 2 
are too blue in i — z to have been sc 
cuts usually applied in optical searches 
quasars. Currently, there is an apparent 

dates at 2 > 6.4 
e radio loudest 
8035 (2 = 5.68) 
78.0 ± 40.5 and 
R = 91.4 + 8.8. 
lft end of the i- 
; 5.7) and they 
lected by color 
or high-redshift 
lack of quasars 

at 5.2 < 2 < 5.7 (see McGreer et al. 

2013 Banados 


between optical colors of quasars and M dwarfs which 


a re th e most numerous stars in the Galaxy (Rojas-Ayala 


et al. 2014). The identification of these two radio-loud 


quasars m our extended selection criteria implies that 
there must be a significant number of radio-quiet quasars 
at these redshifts that are just being missed by standard 
selection criteria. The use of additional wavelength infor¬ 
mation, for example using WISE photometry, might help 
to find quasars in this still unexplored redshift range. 

We inspect all the 1.4 GHz FIRST images of the known 
quasars that are in the FIRST survey footprint but not 
in the catalog. Based on this inspection, we identify two 


additional 2 ~ 6 radio-loud quasars which are detected 
at S/N > 3 and therefore would benefit from deeper radio 
imaging. 

We highlight the importance of infrared photometry 
(e.g., from Spitzer or WISE) for z > 5.5 quasars in or¬ 
der to have an accurate measurement of the rest-frame 
4400 A luminosity which allows us to robustly classify 
quasars as radio-loud or radio-quiet. By using Spitzer 
photometry we reclassify the quasar J0203+0012 at z = 
5.72 as radio-quiet (R = 4.3 ± 0.5). This quasar was 
previously classified as radio-loud by estimating its rest 


fra me 44 00 A lumi nosity from the magnitude at 1450 A 


(Wang et al. 

2008 

). The estimate based on an infrared 

proxy is muc 

h better than the one based on the 777-1450 


proxy because less extrapolation is needed. Thus, the es¬ 
timated rest-frame 4400 A luminosity is less affected by 
spectral energy distribution assumptions. 

We compile all the z > 5.5 quasars having 1.4 GHz 
data in the l itera ture and, by making simple assumptions 
(see Section 6.3), we find that the radio-loud fraction of 
quasars at 2 ~ 6 is between 6 % and 19%. We also es¬ 
timate the radio-loud fraction using the Kaplan-Meier 
estimator which takes into account both radio detections 
and upper limits, obtaining RLF = 8 .ll|' 2 %. This frac¬ 
tion suggests no strong evolution of the radio-loud frac¬ 
tion with redshift. This result contrasts with some lower 
redshift studies that show a decreas e of the ra dio-lou d 
fraction of quasars with redshift (e.g., Jiang et al. 2007). 

The study of the RLF of quasars at z > 6 has the 
potential to give a definitive answer to the issue of a 
possible evolution of the RLF of quasars with redshift. 
For this, a homogeneous radio (and infrared) follow-up 
of a well-defined sample of z ~ 6 quasars (or z ~ 7 
when more of these objects are discovered) selected in 
a consistent manner is crucial to test whether there is 
evolution in the RLF of quasars. 
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Figure 2. Rest-frame 5 GHz radio luminosity versus rest-frame 4400A optical luminosity for all z > 5.5 quasars in the literature observed 
at 1.4 GHz (see Table [ 3 ). Filled symbols are 1.4 GHz detections and arrows represent 3a upper limits in the rest-frame 5GHz luminosity. 
The dashed lines represent the radio-to-optical ratios ( R ) of 100 and 10 (separation of radio-loud and radio-quiet). The red circles are 
the two new Pan-STARRSl radio-loud quasars presented here. These quasars are among the most radio loud quasars at thes e redshifts. 
The red downward triangles are two quasars that we identify as radio-loud by inspecting their FIRST images (see Section |6.2| >. The blue 
downward triangles are upper limits for 28 quasars that are located in the FIRST survey footprint. The blue upward triangle represents 
the location of the stack of the 28 undetected quasars in the FIRST survey assuming their average redshift of 6. The gray cross represents 
the position of J0203+0012 if its T 4400 A calculated from extrapolating 1711450 and it is connected by a dotted line with the position 
determined by estimating T 44 qq^ from the IRAC 3.6/im photometry (see Sec tion | 6.1|l. The horizontal gray line denotes the division of 
radio-loud and radio-quiet quasars employed by |Jiang et al. (200?|l (see Section |G.4p 
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Figure 3. 1.4 GHz postage stamp images from the FIRST Survey (typical resolution of 5") of the quasars J1609+3041 (left) and 

J2053+0047 (right). The blue circles have a radius of 3" and are centered on the quasar optical positions. These positions have J2000 
coordinates corresponding to (R.A., Decl.) = (16:09:37.28, +30:41:47.7) for J1609+3041, and (R.A., Decl.) = (20:53:21.77, +00:47:06.8) 
for J2053+0047. The red crosses show the 1.4 GHz peak pixel. The black solid (dashed) lines are the positive (negative) 1,2,3c -contours. 
J1609+3041 and J2053+0047 have a S/N of 3.5 and 3.0, respectively. Using these radio flux densities, both quasars are classified as 
radio-loud with R = 28.3 ± 8.6 and R = 44.1 ± 18.7, respectively (see Table |3j). 
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Table 1 

Candidates After Visual Inspection 


Candidate 

R.A. 

Decl. 

9P1 

rpi 

2 pi 

zpi 

ypi 

&1.4 GHz,peak 
(mJy) 

Prio. a 


(J2000) 

(J2000) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 


PSO J141.7159+59.5142 

09:26:51.83 

+59:30:51.2 

> 23.69 

> 23.39 

> 23.27 

> 22.71 

20.64 + 0.14 

1.11 ± 0.15 

L 

PSO J172.3556+18.7734 

11:29:25.36 

+18:46:24.3 

> 23.64 

> 23.41 

> 23.26 

> 22.78 

20.68 + 0.16 

1.02 ± 0.14 

L 

PSO J044.9329-02.9977 

02:59:43.91 

-02:59:51.9 

> 22.81 

> 23.33 

22.88 ± 0.28 

21.16 + 0.08 

20.89 + 0.17 

5.86 ± 0.08 

H 

PSO J049.0958-06.8564 

03:16:23.00 

-06:51:23.2 

> 22.57 

> 23.11 

23.07 + 0.28 

21.37 + 0.10 

20.76 + 0.17 

2.77 + 0.14 

H 

PSO J055.4244-00.8035 

03:41:41.86 

-00:48:12.7 

> 22.88 

> 23.08 

> 21.54 

20.28 + 0.05 

20.27 + 0.10 

2.14 + 0.14 

H 

PSO J106.7475+40.4145 

07:06:59.40 

+40:24:52.3 

> 23.57 

> 23.44 

> 22.73 

21.39 + 0.10 

20.85 + 0.17 

1.37 + 0.13 

L 

PSO J114.6345+25.6724 

07:38:32.30 

+25:40:20.8 

> 23.55 

22.95 ± 0.14 

22.48 + 0.25 

20.89 + 0.10 

20.62 + 0.15 

6.75 ± 0.13 

H 

PSO J135.3860+16.2518 

09:01:32.65 

+16:15:06.8 

23.61 ± 0.24 

> 23.97 

22.38 + 0.17 

20.67 + 0.05 

20.82 + 0.15 

3.04 + 0.14 

H 

PSO J164.9800+07.4459 

10:59:55.22 

+07:26:45.5 

> 23.08 

> 22.68 

21.71 ± 0.14 

20.17 + 0.05 

> 21.39 

3.35 ± 0.14 

H 

PSO J208.4897+11.8071 

13:53:57.54 

+11:48:25.6 

> 23.46 

23.27 + 0.28 

22.73 + 0.15 

20.96 + 0.09 

20.79 + 0.14 

2.28 + 0.13 

H 

PSO J238.0370-03.5494 

15:52:08.89 

-03:32:58.0 

> 23.57 

> 23.6 

22.85 + 0.33 

21.15 + 0.07 

20.84 + 0.19 

6.01 ± 0.15 

H 

PSO J354.6110+04.9453 

23:38:26.65 

+04:56:43.3 

> 23.46 

23.01 ± 0.19 

22.77 + 0.18 

20.88 + 0.09 

20.98 + 0.29 

6.44 + 0.13 

H 


Note. — The two entries at the top are 2 -dropouts and the ten at the bottom are i-dropouts. The lower limits correspond to 3cr limiting magnitudes. 
a Priorities. H: High. L: Low 
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Table 2 

High Priority z-dropout Candidates Follow-up. 


Candidate 

9GROND 

(mag) 

J’GROND 

(mag) 

*GR.OND 

(mag) 

^GROND 

(mag) 

Jgrond 

(mag) 

^GROND 

(mag) 

J'GROND 

(mag) 

Jntt 

(mag) 

Zntt 

(mag) 

Jntt 

(mag) 

Note a 

PSO J044.9329-02.9977 

> 23.72 

23.37 + 0.29 

22.20 + 0.18 

22.14 + 0.23 

20.35 + 0.17 

19.92 ± 0.17 

19.20 + 0.28 

— 

- 

- 

3 

PSO J049.0958-06.8564 

24.47 + 0.33 

23.90 ± 0.23 

23.12 ± 0.26 

21.98 + 0.12 

20.52 + 0.15 

20.01 ± 0.18 

19.19 + 0.24 

- 

- 

- 

3 

PSO J055.4244-00.8035 

> 23.73 

> 23.77 

22.16 + 0.18 

20.58 + 0.05 

20.08 + 0.16 

20.03 ± 0.22 

> 19.09 

- 

- 

- 

1 

PSO J114.6345+25.6724 

- 

- 

- 

- 

- 

- 

- 

22.04 + 0.18 

20.95 + 0.28 

20.48 + 0.11 

2 

PSO J135.3860+16.2518 

> 24.57 

24.32 ± 0.36 

22.70 + 0.18 

20.85 + 0.04 

20.30 + 0.12 

20.91 ± 0.33 

> 19.71 

- 

- 

- 

1 

PSO J164.9800+07.4459 

23.66 + 0.15 

22.19 + 0.06 

21.15 ± 0.05 

20.49 + 0.02 

19.96 + 0.10 

19.72 ± 0.13 

19.45 + 0.27 

20.64 + 0.04 

20.27 + 0.04 

- 

3 

PSO J208.4897+11.8071 

> 24.50 

23.21 ± 0.18 

22.08 + 0.25 

21.66 + 0.11 

20.53 + 0.23 

19.83 + 0.19 

> 19.46 

- 

- 

- 

3 

PSO J238.0370-03.5494 

> 24.47 

23.44 + 0.18 

22.03 + 0.11 

21.32 + 0.07 

19.99 + 0.11 

19.74 + 0.16 

> 19.70 

- 

- 

- 

3 

PSO J354.6110+04.9453 











2 


Note. — The magnitude lower limits correspond to 3cr limiting magnitudes. 

a 1: ^ > 5.5 quasar spectroscopically confirmed in this work. 2: Not a z > 5.5 quasar based on follow-up spectroscopy. 3: Not a z > 5.5 quasar based on follow-up photometry. 
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Table 3 

Data and Derived Parameters of the 2 : > 5.5 Quasars with 1.4GHz Data in the Literature. 


Quasar 

z 

Si .4 GHz,peak 

(^Jy) 

Ref.(z,1.4GHz) 

mi450 a 
(mag) 

Wi b 
(mag) 

53.6/im C 

(mag) 

l°g £5 GHz 

( L o) 

lo gf4400A d 

(£©) 

R d 

J0002+2550 

5.82 

89 ± 14 

1,2 

19.0 

18.86 ±0.06 

18.71 ± 0.02 

8.03 ±0.07 

13.05 ±0.01 

1.3 ±0.2 

J0005-0006 

5.85 

< 390 

3,4 

20.2 

20.00 ±0.16 

20.10 ±0.03 

< 8.7 

12.49 ±0.01 

< 20.9 

J0033-0125 

6.13 

< 57 

5,4 

21.8 

20.99 ±0.40 

- 

< 7.9 

12.19 ±0.15 

< 6.8 

J0203+0012 

5.72 

195 ± 22 

6,4 

21.0 

19.39 ±0.09 

19.14 ± 0.03 

8.36 ±0.05 

12.86 ±0.01 

4.3 ±0.5 

J0303-0019 

6.08 

< 186 

7,4 

21.3 

- 

20.24 ± 0.04 

< 8.4 

12.47 ± 0.01 

< 11.4 

PJ055-00 

5.68 

2140 ± 137 

8,9 

20.4 

20.62 ±0.26 

- 

9.39 ±0.03 

12.27 ± 0.09 

178.0 ± 40.5 

J0353+0104 

6.049 

< 57 

10,4 

20.2 

19.34 ± 0.09 

19.44 ± 0.04 

< 7.9 

12.79 ±0.01 

< 1.7 

J0818+1722 

6.02 

123 ± 12 

1,2 

19.3 

- 

18.35 ±0.01 

8.20 ±0.04 

13.22 ±0.01 

1.3 ±0.1 

J0836+0054 

5.81 

1740 ± 40 

3,2 

18.8 

18.00 ±0.04 

17.87 ± 0.01 

9.32 ±0.01 

13.38 ±0.01 

11.9 ±0.3 

J0840+5624 

5.8441 

< 27 

11,2 

20.0 

19.46 ±0.14 

19.53 ±0.02 

< 7.5 

12.72 ±0.01 

< 0.9 

J0841+2905 

5.98 

< 81 

1,4 

19.6 

19.91 ± 0.16 

19.74 ± 0.05 

< 8.0 

12.65 ±0.02 

< 3.1 

J0842+1218 

6.08 

< 57 

12,4 

19.6 

- 

19.13 ±0.01 

< 7.9 

12.91 ± 0.01 

< 1.3 

PJ135+16 

5.63 

3040 ± 145 

8,9 

20.6 

19.51 ± 0.11 

- 

9.53 ±0.02 

12.71 ± 0.04 

91.4 ±8.8 

J0927+2001 

5.7722 

50 ± 11 

13,4 

19.9 

19.40 ±0.11 

19.72 ±0.05 

7.77 ± 0.10 

12.63 ±0.02 

1.9 ±0.4 

J1030+0524 

6.308 

< 60 

3,2 

19.7 

19.28 ±0.09 

19.23 ±0.04 

< 7.9 

12.91 ± 0.02 

< 1.5 

J1044-0125 

5.7847 

< 72 

14,2 

19.2 

19.05 ±0.07 

18.84 ± 0.02 

< 7.9 

12.99 ±0.01 

< 1.2 

J1048+4637 

6.2284 

< 33 

1,2 

19.2 

19.05 ±0.06 

18.80 ± 0.01 

< 7.7 

13.07 ± 0.01 

< 0.5 

J1120+0641 

7.0842 

< 23 

15,16 

20.4 

19.61 ± 0.11 

19.39 ±0.03 

< 7.6 

12.93 ±0.01 

< 0.7 

J1137+3549 

6.03 

< 51 

1,2 

19.6 

19.16 ±0.07 

19.09 ±0.03 

< 7.8 

12.92 ±0.01 

< 1.1 

J1148+5251 

6.4189 

55 ± 12 

1,17 

19.0 

18.67 ± 0.05 

18.57 ± 0.02 

7.91 ± 0.09 

13.18 ±0.01 

0.7 ± 0.2 

J1250+3130 

6.15 

< 63 

1,2 

19.6 

19.11 ± 0.07 

19.09 ±0.01 

< 7.9 

12.94 ± 0.01 

< 1.3 

J1306+0356 

6.016 

< 63 

3,2 

19.6 

19.57 ± 0.10 

19.24 ± 0.04 

< 7.9 

12.86 ±0.02 

< 1.5 

J1319+0950 

6.133 

64 ± 17 

14,18 

19.6 

19.73 ±0.11 

- 

7.94 ± 0.12 

12.69 ±0.04 

2.4 ± 0.7 

J1335+3533 

5.9012 

35 ± 10 

1,2 

19.9 

19.41 ± 0.07 

19.35 ± 0.02 

7.64 ± 0.12 

12.80 ±0.01 

0.9 ±0.3 

J1411+1217 

5.904 

61 ± 16 

3,2 

20.0 

19.29 ±0.07 

19.05 ±0.02 

7.88 ±0.11 

12.92 ±0.01 

1.2 ±0.3 

J1425±3254 e 

5.8918 

< 60 

1,4 

20.6 

19.67 ± 0.08 e 

20.36 ± 0.06 e 

< 7.9 

12.39 ±0.02 

< 4.1 e 

J1427+3312 

6.12 

1730 ± 131 

19,9 

20.3 

19.52 ±0.08 

19.49 ±0.02 

9.37 ± 0.03 

12.77 ±0.01 

53.3 ±4.1 

J1429+5447 

6.1831 

2930 ± 152 

18,9 

20.9 

19.73 ±0.08 

- 

9.60 ±0.02 

12.70 ±0.03 

109.2 ± 8.9 

J1436+5007 

5.85 

< 48 

1,2 

20.2 

19.87 ± 0.09 

19.79 ±0.02 

< 7.8 

12.62 ±0.01 

< 1.9 

J1509-1749 

6.121 

< 54 

20,18 

19.8 

- 

- 

< 7.9 

12.91 ± 0.11 

< 1.2 

J1602+4228 

6.09 

60 ± 15 

1,2 

19.9 

18.75 ±0.04 

18.57 ± 0.02 

7.90 ±0.11 

13.14 ±0.01 

0.8 ±0.2 

J1609+3041 

6.14 

484 ± 137 

21,22 

20.9 

20.22 ±0.14 

- 

8.82 ±0.12 

12.50 ±0.05 

28.3 ±8.6 

J1621+5155 

5.71 

< 63 

4,4 

19.9 

18.35 ±0.03 

- 

< 7.9 

13.18 ±0.01 

< 0.7 

J1623+3112 

6.26 

< 93 

18,2 

20.1 

19.22 ±0.06 

19.23 ±0.03 

< 8.1 

12.90 ±0.01 

< 2.3 

J1630+4012 

6.065 

< 45 

1,4 

20.6 

20.19 ±0.12 

19.98 ±0.06 

< 7.8 

12.57 ± 0.02 

< 2.2 

J1641+3755 

6.047 

< 96 

20,4 

20.6 

- 

- 

< 8.1 

12.58 ±0.12 

< 4.6 

J2053+0047 

5.92 

434 ± 143 

23,22 

21.2 

20.82 ±0.32 

- 

8.74 ± 0.14 

12.23 ±0.12 

44.1 ± 18.7 

J2054-0005 

6.0391 

< 69 

14,4 

20.6 

- 

- 

< 8.0 

12.59 ±0.12 

< 3.2 

J2147+0107 

5.81 

< 54 

23,18 

21.6 

20.33 ±0.20 

- 

< 7.8 

12.41 ± 0.07 

< 3.5 

J2228+0110 

5.95 

310 ± 57 

24,24 

22.2 

- 

- 

8.59 ±0.08 

11.94 ± 0.13 

61.3 ±20.9 

J2307+0031 

5.87 

< 51 

23,18 

21.7 

19.78 ±0.13 

- 

< 7.8 

12.64 ± 0.05 

< 2.0 

J2315-0023 

6.117 

< 48 

10,4 

21.3 

20.26 ±0.20 

20.10 ±0.03 

< 7.8 

12.53 ±0.01 

< 2.6 

J2329-0301 

6.417 

< 66 

20,4 

21.6 



< 8.0 

12.21 ± 0.12 

< 8.3 


References. — (1 ) iCa rilli et al.| 
Mortlock et al. (2009), (7) Kurk etf 
(V2) De Rosa et ai. ( 201ljh~ (13) |UarT. 


( 201( 51) , ( 2) I Wang et al. I d2007h , (3) iK urk et a l. (2007), (4) Wang et al. (200 81, (5) IWjjjgttT 
al. (2009), (8) This Work, (9) FIRST Becker et al. (1995), (10) Jiang et al. (2008), (II) W. 
arflli et ak (2007), (14) Wang et al. (2U13), (15) Vcnemans et al. (2012), (1 b) |lVlomjiarie' 


(12) De Rosa et al. (201 l b , (13 ) (Jarilli ct al. (200 71), (14) I Wang e t ah (2U13), (15) Vcnemans et alf Tl 
Carilli et al. (2004b), (18) Wang ct al. (2011k (19) McGrecr et al. (2006), (20) Willott et al. (2010a), 
Work: new radio identification (see Section |b.2[ ), (23) Jiang ct al. ( 20 q 5]), (24) [Zeimann ct al. ( 2Ull| ) 


et al. 1 (|200T[> ■ (6) 

ling et al!]l 2010 ), 

et al.|(|2U14rpT) 
(21) Warren et al. (in prep.), (22) This 


Note. — Reported upper limits correspond to 3cr. 

a All 777.1450 are taken from|Calura et al.|(|2014|) except for J2228+0110 for which is taken from|Zeimann et al.|(|2011|, for J1609+3041 for which 
is calculated from its ypi band (Banados et al. in prep.), and for PJ055-00 and PJ135+16 for winch are calculated in this work as in|Banados| 
|et al.| ( [2014| . 

b All W\ measurements have S/N > 2.5. Magnitudes are taken from the main ALLWISE source catalog with exception of J0033-0125, PJ055-00, 
and J2053+0047, for which are t aken from the ALLW ISE reject table. 

c Sp.p.jr rr, measurements are f rom|Leipski et al.|(|2014[) with exception of J1120+0641 and J1425+3254 for which are taken from | Bar nett et al.| 
( |2015f ) and |Cool et al.| ( |2006[ ) , respectively. 

d log L 4400 A an d f v 4400 A i n R — fv ,5 GHz/ fv 4400 A are based on S 3 . 6 p, m measurements when available, otherwise from W\. If the quasar does 
not have S3. 6// m nor W\ data, the quantities are extrapolated from 7711450 ( see text in Section |5.2| ). 

e We note a large discrepancy between the reported WISE and Spitzer magnitudes for J1425+5254. If W± is used instead of S 3 . 6 p,m to estimate 
fv, 4400 A> R would be < 2.1. 



















































































